arXiv:1506.02942v2 [physics.plasm-ph] 19Jun2015 


Transport Barrier generation 

at the interface of regions with different zonal flows dynamics 

C. NorsciniQPh.Ghendrih, T. Cartier-Michaud, G. Dif-Pradalier, X.Garbet, N.Nace, Y. Sarazin, and P. Tamain 

CEA, IRFM, F-13108 St. Paul-lez-Durance cedex, France 
(Dated: June 22, 2015) 

A novel and generic understanding of spontaneous generation of transport barriers and zonation 
regimes in turbulent self-organization is presented. It associates the barrier onset to the development 
of a spectral gap between large scale flows and turbulence modes leading to a zonation regime. A 
robust barrier builds-up at the interface of such a region and a neighboring one with reduced zonal 
flow generation. This more complex and generic transition paradigm could fit the numerous and 
sometimes conflicting observations as in fusion plasma experiments. Barrier relaxation by bursts 
of turbulence regenerate the zonal flows that are eroded by viscous (collisional-like) damping. The 
duration of the quiescent phase between the quasi-periodic relaxation events is governed by this 
damping process, hence the barrier collision frequency for fusion plasmas. 
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Turbulence self-organization plays a major role in 
transport properties within stratified media in geo¬ 
physics [T], astrophysics [2] as well as laboratory plasma 
dedicated to magnetic fusion research 13 a. An 
outstanding mechanism is the transition from a fully 
turbulent to the so-called zonation regime [5], where 
large scale anisotropic flows, the zonal flows, appear to 
undergo condensation to a regular pattern [6j |7] . These 
regimes are investigated via reduced 2D models of the 
Hagesawa-Mima type with prescribed energy injection 
in the spectra of otherwise decaying turbulence. They 
are regarded as a possible mechanism for the generation 
and maintenance of stable systems of alternating zonal 
jets as observed in the deep terrestrial oceans and in the 
atmosphere of giant planets [8]. In the fusion plasma 
framework, self-organized zonal flows are considered to 
govern transport barriers [3 HO]. To study such relation, 
particle (or heat) transport has to be taken in account 
in the 2D model. 

Without loss of generality, we address such problem 
via a model comparable to the Rayleigh-Benard (RB) 
instability, where self-consistent turbulence drives both 
the transport properties and the energy injection into 
the turbulent modes via the interchange (buoyancy) 
instability. We report here the formation of a transport 
barrier, or pedestal, localized at the interface between 
two regions with different zonal flow damping capability. 
For magnetic fusion plasma, that correspond to the edge 
region, where the magnetic surfaces are closed, and the 
Scrape-Off Layer (SOL) region [11], where the magnetic 
surfaces intersect wall components. This interface is 
referred to as the separatrix. Of interest here is the 
role of large scale damping in controlling the transport 
regime [laini and consequently the transport barrier. 

The 2D fluid model of the interchange instability 
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m m, is similar to several models addressing turbu¬ 
lence generic properties. It is simplified to only retain 
the key features required for the barrier generation. It 
stands for particle and charge conservation equations in 
an isothermal plasma in the cold ion limit: 

dtn + [0, n] — D±A±n = S — T (la) 

dtW + [0, W] - AxIF + gdyu = J (lb) 

where n is the particle density and W the vorticity, 
W = Ax4>. Time and space are normalized by the in¬ 
verse of the cyclotron frequency and the Larmor radius 
respectively. Space is reduced to 2D transverse to the 
magnetic field assuming symmetry along the magnetic 
field. The y coordinate is an angle, the poloidal angle 
for magnetic confinement, the x coordinate in the ra¬ 
dial direction that extends from the source S', localized 
in the edge, to the sink T in the SOL. The separatrix is 
defined at a given xsep position. Constraints due to the 
physics along the magnetic field govern volumetric source 
terms: the parallel divergence of the electric current J, 
the particle sink T specific of the SOL region. Convective 
turbulent transport [4>,/] = 5^(/(-5^4>)) Fdy{f{dx^)) 
competes with weak diffusive transport with coefficient 
D for the particles and u for the vorticity. The buoyancy 
force governs the plasma interchange instability and 
thus the self consistently energy injection into the tur¬ 
bulent spectrum. 0 is then the electric potential and W 
is related to the polarization charge. For standard fluids 
(/) is the current function. In the RB model n stands for 
the temperature variation. The change in field line prop¬ 
erties at the separatrix is taken into account by a mask 
function x(x), x(x > xsep) = 0 and x(x < xsep) = 1 
such that J = — x(^) < ^ >y),where is the 

normalized conductivity and < / >^ is the ^—average 
of /, similarly F = {1 — x{x))anU where dn is the par¬ 
ticle lifetime in the SOL. The edge and SOL difference 
is twofold. First the particle loss is localized in the SOL 
region. This has rather little effect on the turbulence 
but organizes the overall stratification of the system in 
the X direction. Second, the change in J modifies the 
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evolution equation for the zonal flow Vz = dx < (j) >y 
The latter is governed by a balance between the non 
linear Reynolds stress source, and the loss terms: vis¬ 
cous damping at small scales, sink term < J >y at large 
scales. The edge constraint < J >y= 0 favors large scale 
zonal flow structures. Conversely, in the SOL, the current 
loss < J >y= (J(j) < here linearized for simplicity, 

damps these zonal flow structures. We investigate near 
marginal cases for a prescribed density gradient length 
1/Ln- The linear growth rate is given by: 
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{VCF^) 


1/2 


Lr ky 1 




LiK) \ Ln 5 , 


( 2 a) 


iii= — ; L{K) = {l + llS,)K‘^ + l/K‘^ ( 2 b) 

V(T(f) 


where K = k/k^ k"^ = (Tc^jv and Sc = yjD is the 
Schmidt number. This expression yields three key 
aspects of the near marginal regime, the order of 
magnitude of 7 ^, the threshold depending 

on the Rhine scale Lr [SJ [ 12 ] such that LrIL^ must 
exceed a given function of kx and ky^ a localization 
function L{K) that favors K = {1 as most 

unstable mode. The latter effect is governed by the 
balance between the homogeneous damping rates, at 
small scale due to diffusion at large scale due 

to parallel currents ~ The threshold favors the 

modes with small values of kx^ excluding kx = 0 to en¬ 
sure momentum conservation by the zonal flow structure. 


SOL-transport with the present model is comparable 
to that previously reported when only the SOL region 
was addressed [HUS]. It is characterized by avalanche 
transport, hence ballistic propagation of fronts and holes 
im US], the so-called ’blobs’ routinely observed in ex¬ 
periments. Density and potential fluctuations are large 
[TTJ [T5| but the mean value of the latter weakly departs 
from the equilibrium value, namely the plasma float¬ 
ing potential set at zero for convenience in this isother¬ 
mal case. Conversely, transport in the edge region ap¬ 
pears to be controlled by zonal flows. These generate 
transport barriers (TBs), where the transverse turbu¬ 
lent avalanches are damped and therefore where diffusive 
transport governs a larger fraction of the particle out- 
flux. At the edge and SOL interface a barrier is readily 
observed, flgjTJQ. Three different average density pro¬ 
files are compared on figjTJQ where only the position of 
the separatrix xsep is modified in the simulations, i.e. 
xsep = 0 -b Xa (continuous line), xsep = 0-8 (dashed 
line) and xsep = dash-dot line). One readily observes 
3 key features: the density e-folding length in the SOL 
is unchanged, the pedestal region with increased density 
gradient at the separatrix extends both in the edge and 
in the SOL region, other regions in the density profile 
exhibits corrugations (enhanced gradients) which drive 
staircase like profiles IS US). 

In order to quantify the extent of the barrier in the 
two regions, we shift the profiles to the same separa- 



FIG. I: Q Density pedestal for different positions of 
the separatrix (see text), 0 density gradient profiles 
with radial shift, Q zoom on the pedestal region 


trix position and compare the density gradient profiles 
1 / L/fi — SxT^ / ri, The strong density drop ob¬ 

served in the pedestal leads to a marked peak in the l/L^ 
profile localized at the separatrix. The barrier extends 
into both the edge and SOL region and its width is ob¬ 
served to range between 5 % and 10 % of the minor radius 
figjTJQ. Regions with large zonal flows shear are corre¬ 
lated with the corrugation of the profile, see figj2]Q&Q. 
They are characterized by a stopping capability of most 
of the avalanches both overdense from uphill and holes 
from downhill. 

The simulation is characterized by a slow reorganization 
of the zonal flow pattern as readily observed on the con¬ 
tour plot of the zonal flow shear superimposed on the 2D 
plot of 1/Lnj fig(2]0. 

Two features are outstanding, the evolution towards a 
dipolar structure of < ^ in the edge region while the 
structure of the maximum value at the separatrix weakly 
evolves. When a statistical steady state is reached, one 
can average the profiles over time, figllE and figHii- 
One can see that the total flux T^ot is radially constant in 
the edge and decays in the SOL. The turbulent contribu¬ 
tion Tfurb exhibits well defined minima, and conversely 
large values of the diffusive flux Tdif^ of l/Ln and of 
the zonal flow shear. Narrow regions with strong tur¬ 
bulent transport are localized in the vicinity of the zero 
zonal flow shear layers. To identify the transport barri¬ 
ers, one defines the ratio between the y-averaged particle 
fluxes = Tturb/^tot [IZ]- Rb varics between 0 and I 
in steady state and is a measure of the effectiveness of 
the barrier in reducing turbulent transport. One readily 
observes on figHE) that Rb changes in time and space (x- 
direction). In space, one recovers the dipolar structure 
with four transport barrier regions, a pedestal at the sep- 
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FIG. 2: ^ E X B shear (red contour) and l/L^ (gray 
colored) in function of time and radius,® the time and 
y-averaged shear ( black line) and 1/Ln ^ed circles), Q 
parameter used to determine the transport barriers, 
varies between 1 (white) and 0 (black), total flux 
(black line), turbulent flux (blue diamonds) and 
diffusive flux (red circles) in the radial direction. 


aratrix that is relatively narrow, two broad barriers in 
the edge and finally a small transport barrier towards 
the source region that is strongly linked to the boundary 
conditions of the model. In time, one can observe quasi- 
periodic relaxation events characterized by strong turbu¬ 
lent transport across all the barriers. While these events 
are globally quasi-periodic, the detailed time evolution 
is specific of each event made of consecutive avalanches 
that do not extend throughout the edge region [THl HO]. 
These events are also correlated with the large transport 
bursts in the SOL region, see figHE- As can be observed 
on the time traces, the edge region exhibits a sawtooth 
structure corresponding to storage by the TB while the 
SOL region exhibits a pulse-like variation since the SOL 
acts as the sink for the particles released at each relax¬ 
ation, figj3]|^&0. 

The generation of a barrier, associated in this model to 
a corrugation structure in the edge region, is observed for 
values of the control parameters that are close to the in¬ 
terchange marginal values. The interplay zonal-turbulent 
modes can be addressed in the framework of nonlinear 
three mode coupling, (psi^^ky), (j)t{i^^—ky)^ re¬ 

spectively the zonal flow, a streamer and a more homo¬ 
geneous turbulent mode [iniiin]- Setting the interchange 
instability to zero, ^ = 0 , we address the dispersion rela¬ 
tion in two cases (1) Z-flow generation with growth rate 
7 i, finite streamer (j)s mode as reference equilibrium and 
two coupled modes as perturbations (j)z^ conversely, 
( 2 ) zonal flow saturation with growth rate 72 , finite zonal 
equilibrium mode (j)z^ and perturbations 0^, (j)t. The dis- 




FIG. 3: Q l/Ln evolution in time and radial direction: 
the different transport reorganization in SOL and edge 
region, evolution of 1/Ln in time at given radial 
positions in the edge (|^, and SOL 0 


persion relations for case ( 1 ) and ( 2 ) are: 

(71 + lz){li + It) = (3a) 

(72 +7t)(72 +7s) = (3b) 

where ^z — Is — J^ky + cr/ky and jt = E o /k^ ^ 

k^ — ^ k^. The coupling term is Vi = {k‘^ ~ 

It governs a symmetric necessary condition for positive 
growth rates is obtained: in case (a) ky > k,‘^; in case (b) 

kl < 

Gonsidering the turbulence spectrum in ky^ we define 
three regions: the zonal flow (Z) ky = 0, the Big (B) 
\ky\ < |/i:| and Small (S) turbulent structures \ky\ > 
Region S is the source of zonal flows while B acts as a 
sink and energy loss via the large scale damping term 
< 70 . The viscosity term is a generic version of large scale 
damping, including collisional damping, provided 7 ^ 0 , 
that actually controls the linear damping of the zonal 
flows [ 21 ]. The existence of the B-mode region is es¬ 
sential for the transition from a turbulent regime to the 
zonation regime characterized by the development of a 
non homogeneous spectrum, \ky\ < \kx\ ^ hz. Further¬ 
more, the shearing capability of the k^ = tz zonal flow is 
small for the B-modes since \ky\ < For ^ 7 ^ 0, the 
energy injection in the spectrum is governed by the in¬ 
terchange dX ky '^k and the inverse cascade is controlled 
by the Rhine scale eq. (§• The development of the 
B-mode gap between turbulence and zonal flows is thus 
constrained by tZR > tz^ where hZR = min(^, I/Lr). 

The spectrum of the electric potential in simulations 
in the zonation regime, figj40, is characterized by a gap 
between S and Z modes, hence with weak B-mode tur¬ 
bulence. The S-region transfers energy via non linear 
coupling towards the Z-flows, which tends to store the 
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FIG. 4: Q Density gradient profile and 0 electric 
potential spectrum : zonation regime (black line) and 
turbulent regime (red circles) 


energy and quench the turbulent transport dominated 
by the S-modes. The Z-mode is gradually damped by 
the viscosity damping while the TBs gradients build-up 
until a new relaxation event is triggered during which 
the gradients relax and Z-flows are regenerated by the 
S-mode turbulent activity. When the gap between S and 
Z regions is reduced, hence increasing the B-mode am¬ 
plitude, the B-mode turbulent activity is not affected by 
Z-ffow shearing and the Z-mode energy can be transferred 
back to the turbulent modes, both the B and S modes. 
There is then an increase of the frequency of the turbu¬ 
lent bursts to the point where they cannot be isolated 
from the steady state transport activity and the pedestal 
is smeared-out, figJdjQ. The fluctuations spectra of these 
two regimes present a strong similarity with the exper¬ 
imental observation achieved during this L-H transition 
on the stellarator H-1 [22]. 

The dynamics of this transition between pedestal (High 
confinement) and no-pedestal (Low-confinement) behav¬ 
ior, is captured by the following 0-D predator-prey model 

m- 


dtVn 

Vn 


P 

Vn 


(T + T*) 



(4a) 

(4b) 


^ = 7 .(Vn - Vn*)(l - a,S) - (4c) 


dtB 

B 


= 76 (Vn - Vn*)(l - a^^B) 


(4d) 


The model couples the gradients Vn and the Z, S and 
B modes. The gradient Vn is governed by a balance 
between the source P and transport, both turbulent 
T = B P S and collisional T*. The growth rates for 
S and B modes are 7 ^ and 75 respectively and exhibit a 
threshold in the gradient, Vn*. Non-linear saturation of 
these modes are used, proportional to ag and a^. The 
control parameter of Z-ffow generation by S is p. The 
B-modes act as saturation term on this energy exchange 
via the 1/T dependence. The Z-ffow sink is proportional 
to the viscosity n. The limit cycles of the Z-T interplay 
in simulations and 0-D model are compared in the Z — T 
plane, figjsjQ&Q. Z and T are readily determined by 


eq.Q for the 0-D model. For the simulation output of 
the interchange turbulence, we define Zi = FT(Vz) and 
Tf = FT{R), where FT is the 2D Fourier mode energy 
of the zonal flow velocity W and Reynolds stress R re¬ 
spectively. 



FIG. 5: Limit cycles: Q 2 D simulation, (|^ OD model 

The different positions of the cycles in the Z-T 
plane is determined by the control parameters: from 
a reference case, blue trace, increasing the curvature 
term red trace or decreasing a black trace, figjsjQ. 
Increasing g in the simulation leads to an increase of 
both turbulence and zonal flows, the latter being more 
important, increase of Zi/Ti. A comparable behavior 
is obtained by increasing 7 ^ and 75 , namely the growth 
rate of the interchange instability -governed by g- 
figHB- Decreasing a governs a decrease of both k and 
I/Lr so that the spectrum maximum shifts towards the 
low ky values, reducing the Z-S gap. The turbulence 
amplitude is increased as well as the ratio between B 
and S modes. Consistently, this behavior is recovered 
in the 0-D model by reducing ag and the non-linear 
turbulence saturation, as well as the critical gradient 
Vn*. 

Barrier generation at the interface between regions with 
different zonal flow damping is governed by a zonation 
regime in the plasma edge region. It is triggered by 
a gap in the turbulent spectrum between the wave 
vector of the maximum and zero, the zonal flow wave 
vector. This gap is controlled by the energy injection 
wave vector of the interchange instability and the Rhine 
scale that bounds the inverse cascade. Increasing the 
magnitude of the turbulence drive at given gap reinforces 
the barrier. In the pedestal High confinement regime, 
relaxation bursts of turbulence regenerate the zonal 
flows that are eroded by collisional damping. For fusion 
plasmas the duration of the quiescent phase between the 
quasi-periodic relaxation events is thus governed by the 
ion pedestal collision frequency. 
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